Abstract 



Defocused speckle photography has long been used to measure rota- 
tions of rough surfaces. This paper explains how, by adding a suitably 
perforated mask, some measurement properties, such as range or lateral 
resolution, may be changed at wish. Particularly, the maximum measur- 
able tilt can be significantly raised, although at the expense of poorer 
lateral resolution. Advantages over previously described techniques in- 
clude independent tuning of speckle size and optical system aperture and 
more extended adaptability to different measuring needs. The benefits and 
disadvantages of the new and old techniques are thoroughly compared. 
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1 Introduction 

When a rough surface is illuminated by a laser beam, a granular effect is ob- 
served. This grainy pattern is known as speckle. Though this effect used to 
be considered a nuisance in early coherent photography, it has proved to be a 
valuable tool for the observation of mechanical transformations in solid bodies. 
Different optical systems can be used to obtain a record of speckles. Although 
surface roughness causes their formation, their statistics are generally quite in- 
dependent of the particular microtopography of the surface and are more related 
to the optical set-up for illumination and recording.^ Speckle patterns do not 
offer great information about the rough surface that generates it. However, the 
relation between the speckle patterns taken from a surface before and after a 
transformation may be used to obtain information about it. 

Assuming expanded beam illumination, object in-plane displacements (move- 
ments parallel to the surface under observation) may cause speckle pattern 
translations when the optical system is focused on the object surface, whereas 
out-of-plane rotations (tilts around in-plane axis) will not significantly affect 
the image. On the other hand, if the speckle pattern is recorded at the Fourier 
plane of the object surface, in-plane displacements will not appreciably change 
the speckles, although out-of-plane rotations will displace them.^ When the 
optics is not focused on the object surface, the speckle pattern is said to be 
defocused. Defocused speckle photography is being used today as a mean to 
monitor and measure out-of-plane rotations in rigid solids or their distributions 
in the boundaries of elastic bodies. This paper develops an analysis of measur- 
ing characteristics such as lateral resolution or maximum measurable rotation 
in existing technologies, then presents a new enhancement for them and finally 
compares the measuring power in both systems. 
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Figure 1: Experimental set-up using normal illumination and observation. 

2 Antecedents 

The first paper describing a defocused two-exposure method to measure out- 
of-plane rotations was due to Tiziani,^ and was later extended for vibration 
analysis.'* The principle is simple and can be understood from the following 
consideration: when an initially vertical flat mirror undergoes an out-of-plane 
rotation around its horizontal direction, objects reflected by the mirror seem 
to be displaced upwards (or downwards) proportionally to their distance to 
the mirror. Considering the KirchhofF approximation for the scattering from 
a rough body,^'^ the surface can be thought of as a big set of small mirrors, 
slightly rotated from the mean plane. When the object is tilted, all those small 
mirrors experiment the same out-of-plane rotation and when focusing on a plane 
placed at a distance ( from the surface, the amplitude distributions of the waves 
originated from all those micro-mirrors arc all displaced by the same distance 
which is proportional to The speckle pattern is thus displaced by an amount 
that depends on the tilt and can be used to measure it. Tiziani used coUimated 
illumination and recorded the image at the back focal plane of the lens. 

If normal illumination and observation are used (Fig. 1), the speckle shift is 
given by^ 

= 2//? (1) 

dy = -2fa (2) 

where a, f3 are the (small) rotation angles around the x, y axis of a cartesian 
system placed on the mean plane of the object surface and / is the focal length 
of the recording system. Lateral displacements do not appreciably alter these 
values. 
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CoUimatcd illumination is difficult on extended objects. Gregory considered 
a divergent illumination and showed that when the optical system is focused 
on the plane than contains the image of the point source considering the object 
surface as a mirror, then the speckle shift docs not depend on in-plane dis- 
placements, but only on out-of-plane tilts. ^"^^ Chiang and Juang described a 
method to measure the change in slope by defocused systems. A great number 
of later papers document the use of defocused speckle photography to measure 
in-plane and out-of-plane rotations and strains. This paper describes a new 
methodology that can be applied to almost all of them bringing advantages and 
disadvantages that will be discussed. 

After the double specklegram is obtained by addition (double exposure on 
the same photographic film) or subtraction (using a digital electronic camera), 
a pointwise analysis or a whole field filtering can be performed in order to know 
the lateral distribution of the displacement. ^^'^ The first technique, illustrated 
in Fig. 2(a), uses a narrow laser beam that is made to go through the speck- 
legram. The diffraction pattern is modulated by fringes whose directions are 
perpendicular to the speckle shift and whose spacing is 



where ^ is the distance from the specklegram to the observation plane. 

Whole field filtering, depicted in Fig. 2(b), uses a wide laser beam that goes 
through the specklegram and is filtered by a hole placed in a plane standing 
before a focusing lens that reconstructs an image of the specklegram. Only 
those areas whose modulated diffraction patterns do not have a dark fringe over 
the hole, will appear illuminated on the image plane. Thus, dark fringes in the 
detector can be mapped onto speckle shift values. 

Today CCD cameras replace photographic films and processing of speckle- 
grams is progressively performed through digital techniques,^'' although they 
usually simulate the analogic ones. 

3 Direct illumination system analysis 

We will consider the Tiziani optical system adapted for normal illumination 
represented in Fig.l. A laser source LS emits a light beam, that is expanded 
by BE and spatially filtered by SF. Half of the intensity goes through the beam 
splitter BS and reaches the rough surface S. Part of the light scattered by S is 
reflected by BS and recorded at the back focal plane R of the lens LI. We will 
refer to this set-up as DDSP (direct defocused speckle photography) set-up. The 
surface under observation is supposed to be a rough L x L square, placed at 
a distance d from the objective lens of the optical system, measured along the 
optical axis (dashed in Fig. 1). The aperture diameter of the system is D, its 
focal length is / and the sensing area is a 6 x 6 square whose resolution width 



a = X 




(3) 



is S. 
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Figure 2: (a) Fringes arising from speckle shift around one point, (b) Fringes 
arising from whole field filtering. 



For a defocused recording system, the speckle size is given by^* 

. = i^ (4) 

When a part of the object surface undergoes an out-of-plane rotation 7, the 
light scattered from it experiments a rotation 27 and completely falls off the 
aperture of the system when 

27 > ^ (5) 
Thus the maximum measurable rotation is 



^ (6) 



2dF 
being 

F=^ (7) 

the f-number of the lens. Throughout this paper high apertures will mean low 
f-numbers and vice versa. In order to use the whole recording area, we should 
ensure that 

d--f 
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and accordingly 

'"=21? 

Neglecting diffraction effects (because we are using a dcfocuscd system), the 
diameter of the zone whose scattered light is incident on the same point at the 
detector plane is D, so the lateral resolution A of the system is equal to D 

The relation 

r 6 

A = 27Z (10) 

does not depend on the aperture. 

By properly choosing the system aperture, the technique can improve the 
measuring range or the lateral resolution. Although, according to the equations, 
this versatility may seem feasible, speckle size considerations shrink the choice 
possibilities. Speckle size can not be smaller than the detector cell, because 
averaging destroys the speckle. On the other hand, there should be a large 
number of speckles, in order to accurately track their displacements, what sets 
an upper bound on speckle size. This limits the range of practical values for the 
aperture of the system. Often speckle size is made approximately equal to the 
detector cell size, so that 

and therefore 

1-22/A 

A = (13) 

The main problem with the previously mentioned technique is its short range 
of measurements. The optical aperture of the system must remain small so that 
the speckle size is not smaller than the resolution of the CCD array. This fact 
limits the out-of-plane rotations whose speckle shifts can be recorded by the 
optical system. 



4 Multi-windowed system analysis 

The DDSP set-up described in the previous section can be modified by the 
addition of a suitably perforated mask that enhances the adaptability of the 
system to a wider range of measuring needs. Fig. 3 represents an optical 
system where the laser source LS generates a beam which is expanded by a beam 
expander BE and spatially filtered through the pinhole SF. Then the collimated 
beam is split into multiple collimated beams by a multi-windowcid mask MWM, 
which is shown in Fig. 4(a), with a rectangular array of perforated circular 
holes of radius a {\ <^ a). The beam splitter BS lets half of the light arrive at 
the object surface S and reflects light from S into the aperture of the lens LI. 
Finally a speckle pattern, as represented in Fig. 4(c), is recorded at the back 
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Figure 3: Multi-windowed normal illumination and observation set-up. 



focal plane R of the lens. The focusing and recording optics are implemented 
by an electronic CCD camera. Wc will refer to this set-up as MWDSP (multi- 
windowed defocused speckle photography) set-up. Other arrangements of holes 
are possible, although we will focus on the rectangular array depicted in Fig. 
4(a), which is especially well suited for obtainig a square matrix yij = 7(i A, jA) 
of tilts, where A is the cartesian discretization step of the system. Equilateral 
triangular arrangements, as seen in Fig. 4(b), can also be used and are better 
than rectangular arrays when other parameters such as intensity efficiency or 
lateral resolution are concerned (this is because the non-perforated area and 
the distances between the centers of neighboring holes are smaller). Only the 
rectangular system will be considered here, because for other arrangements the 
analysis is practically identical. 

We will now establish the main relations between relevant parameters of 
the system. In short, the number of illuminated spots on the object surface 
determines the lateral resolution of the system and the width of the detector 
area available for each, which then sets the maximum non-overlapping speckle 
displacement or, equivalently, the tilt range. 

If the object surface under observation and the sensing areas are L x L and 
bxb squares respectively, then, as in DDSP, the distance from the object to the 
camera, measured along the optical axis (dashed in Fig. 3), should be 



In MWDSP, the average speckle size s is given by^^ 

_ 1.22A/ 
* ~ 2a 



(14) 



(15) 
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where a is the narrow beam radius, A is the laser wavelength and / is the focal 
distance of the camera lens. 

Speckle size should be at least as wide as the detector cell, so that 

d<s^a< i^^^ (16) 

When the surface undergoes a tilt of angle 7, speckles are displaced by a 
distance £j 

I = 2^f (17) 

If we consider the maximum tilt to be F, then the maximum displacement will 
be 

Im = 2r/ (18) 

Provided that the radius of the speckle circle is maximum, being N the number 
of rows, we have 

iV<j^ (19) 



The speckle circle radius g is 



where F is the f-number of the camera lens. The non-overlapping condition for 
the recorded speckle circles is given by 

so that 

F>1^ = JL- (22) 
and the usage of the detector area is maximized when 

b 



2Lr 



(23) 



We will assume this optimal aperture to be adopted throughout the rest of this 
paper. 

The non-overlapping condition for the holes in the mask is 

2aN < L (24) 

which leads to 

1.22iVA/ 1.22A/ 

^ = ^ (25) 

or 

s > 1.22XF (26) 

so that speckles in MWDSP are never smaller than in DDSP, provided that 
both systems have the same aperture. 
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(b) 





(c) 



(d) 



Figure 4: (a) and (b): grids of windows that let the laser light reach the ob- 
ject surface in (a) the rectangular and (b) the equilateral symmetries, (c) and 
(d): Speckle patterns at the CCD detector for (c) the rectangular and (d) the 
equilateral arrangements. 
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Table 1: Results for DDSP and MWDSP 



Parameter DDSP MWDSP 

i.m/ 



Speckle size s 1.22AF 
Range F 
Lateral resolution A 
Quotient 



2o 

b b 



2LF 2fN 

L L 
F N 
T b b 



_A 2fL 2,fL 



5 Comparison 

We are now going to compare the measuring range and lateral resolution of the 
DDSP and MWDSP systems. Table 1 summarises the results derived in the 
previous sections. 

In DDSP there is the possibility to change the f-number F of the system. 
As for most popular cameras, 1 < F < 20, then we have 

Fddsp < ^ (27) 

and 

0.05/ < Addsp (28) 

so that the measuring range is only limited by the -| relation. 
In MWDSP, the aperture is related to F through the equation 

2LTmwdsp ^ ^ 

and to the lateral resolution through 

/ 

^MWDSP = (30) 

so that the system aperture relates identically to the measuring parameters in 
DDSP and MWDSP. Thus ^ 

^MWDSP < ^ (31) 

0.05/ < Amwdsp (32) 
Finally we arrive at the relation 

r b 



A 2fL 



(33) 



10 



which mutually limits the measuring range and lateral resolution of the system 
in both DDSP and MWDSP. 

When considering the illumination power, it is straightforward to assert that 

if the aperture is the same in both systems, then there is more efficiency in DDSP 
than in MWDSP. In general, the relation of power at the detector is 

Imwdsp _ nN'^a'^ ^^^^ 
Iddsp 

where Iddsp and Imwdsp are the mean intensities in the DDSP and MWDSP 
systems. Substituting according to the contents of Table 1, yields 

/ 1.22A/ A \^ 

Imwdsp Tra^ \2smwdsp^'^'^^pJ 

= '^-7^^ 72 (35) 



Iddsp 



and simplifying we get 

Imwdsp _ ^ [ sppsp ^ddsp \ ^ 
Iddsp 4 \ smwdsp ^mwdsp J 

or 



Imwdsp _ sddsp Fmwdsp \ 
Iddsp 4 \ smwdsp Fddsp J 



2 



(36) 
(37) 



which, on account of Eq.(26), gives 



Imwdsp < tt ^^^^ 



Iddsp 4 

Finally, substituting for the apertures in Eq.(37) yields 

Imwdsp _^ f sppspTppsp V 
Iddsp 4 \smwdsp^mwdsp J 



(39) 



When both apertures and speckle sizes are equal, the MWDSP power efficiency 
is about 78.5% that of DDSP. If equilateral triangular arrays arc used, this result 
is slightly improved. However, under these conditions, the mcian intensity within 
the speckle circles in the MWDSP system and Idsp are equal. In MWDSP we 
can increase the power at the detector by reducing either the measuring range 
(and accordingly the aperture and lateral resolution) or the speckle size or both, 
in order to get a similar intensity. Nevertheless, when the speckles in MWDSP 
are larger than in DDSP, the intensity in the latter is always greater than in 
the former. This is the main disadvantage of MWDSP that, in our opinion, is 
outweighed by the considerations that follow. 

If we now take into account that the resolution of the detector is 6 there are 
more restraints to be fulfilled. Speckles need to be at least as wide as S so that 
there is no averaging of speckles in the recorded picture. This poses a serious 
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limitation on DDSP since both its speckle size and measurement properties 
depend on F; however, in MWDSP, speckle size only depends on a, so that it 
does not alter the choice range for F and keeps the flexibility of the system 
untouched. Moreover, there is no problem in reducing the value of a (as long as 
we are far from A) making speckles as big as desired. In fact the only restriction 
on MWDSP is that 

2Na <L^ SMWDSP > 1.22XF (40) 

so that the a parameter can always be adjusted to generate speckles whose 
size is 24. 4A which is far from making speckles too large, since generally 6 is 
not much smaller than 24. 4A. In other words, for each value of F there is 
a minimum size of speckle (1.22Ai^) , which is always attainable by choosing 
the maximum admissible a. Even for the maximum value of F, the minimum 
obtainable speckle size is not too large for the usual resolutions of CCD cameras. 
This consideration also applies for DDSP, since 

SDDSP = l-22Ai^ (41) 

is in the _F = 20 limit equal to 24. 4A which is not much bigger than 6. 

Consequently, in practice, speckle size considerations only restraint the mea- 
suring range of DDSP, leaving the lateral resolution limit around 0.05/. Besides, 
the speckle size in MWDSP can be freely raised from the DDSP value to several 
times 5 by reducing the a parameter, without changing the system aperture. 
This possibility makes that the DDSP limitation on the measuring range does 
not hold for MWDSP. 

Table 2 lists the measuring characteristics for both systems when 5 = s. 



Table 2: Results for DDSP (s = 6) and MWDSP 



Parameter DDSP MWDSP 



Speckle size s 6 
Range F 
Lateral resolution A 
Quotient 



2a 

1.22&A b 



2L5 2fN 

1.22/A L 

S N 

T b b 



A 2,fL 2,fL 



Typical values for i5 are in the vicinity of 10-20 microns, so that if L/bK. 5, 
^DDSP is in the order of 10~^ radians for 5 = 2Q microns, whereas Tmwdsp 
can be up to b/{2L) « 10~^ radians. Lateral resolution can be reduced up to 
0.05/. Using a 20mm lens, the minimum A is around 1mm for DDSP and 
MWDSP. 
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6 Conclusion 



By suitably varying the optical aperture, different measuring ranges and lateral 
resolutions are attainable in the DDSP-system, almost exactly as by varying 
the parameter in the MWDSP-system. However, the optical aperture in the 
first system has severe limitations on account of speckle size considerations. On 
the other hand, in the MWDSP-system, speckle size is set by properly choosing 
the a parameter independently from the selection of A^. The main drawback 
in MWDSP is its slightly poorer intensity efficiency, which is outweighed by its 
benefits in those applications where the DDSP measuring range falls too short. 
MWDSP is especially well suited for use with bending plates, where large tilt 
ranges are needed. It could be mentioned, as yet another advantage of the 
MWDSP-system, its adaptive potential. Work on masks with non-equal spacing 
of holes is currently in progress, so that the measuring range and the lateral 
resolution can be changed across the object, achieving a fine lateral resolution 
in some points of limited slope change values and a coarser lateral resolution in 
points with greater expected slope changes. A second enhancement possibility is 
the utilisation of multiple masks (before and after the transformation) in order 
to increase lateral resolution and overcome the limiting relation (see Eq.(33)) 
that holds in both systems. 
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